Surface-sampled simulations of turbulent flow at high Reynolds number

A new approach to turbulence simulation, based on a combination of large-eddy simulation (LES) for the whole flow and an array of non-space-filling quasi-direct numerical simulations (QDNS), which sample the response of near-wall turbulence to large-scale forcing, is proposed and evaluated. The technique overcomes some of the cost limitations of turbulence simulation, since the main flow is treated with a coarse-grid LES, with the equivalent of wall functions supplied by the near-wall sampled QDNS. Two cases are tested, at friction Reynolds number Re$_\tau$=4200 and 20,000. The total grid node count for the first case is less than half a million and less than two million for the second case, with the calculations only requiring a desktop computer. A good agreement with published DNS is found at Re$_\tau$=4200, both in terms of the mean velocity profile and the streamwise velocity fluctuation statistics, which correctly show a substantial increase in near-wall turbulence levels due to a modulation of near-wall streaks by large-scale structures. The trend continues at Re$_\tau$=20,000, in agreement with experiment, which represents one of the major achievements of the new approach. A number of detailed aspects of the model, including numerical resolution, LES-QDNS coupling strategy and sub-grid model are explored. A low level of grid sensitivity is demonstrated for both the QDNS and LES aspects. Since the method does not assume a law of the wall, it can in principle be applied to flows that are out of equilibrium.Comment: Author accepted version. Accepted for publication in the International Journal for Numerical Methods in Fluids on 26 April 201

Stellar Winds on the Main-Sequence II: the Evolution of Rotation and Winds

Aims: We study the evolution of stellar rotation and wind properties for low-mass main-sequence stars. Our aim is to use rotational evolution models to constrain the mass loss rates in stellar winds and to predict how their properties evolve with time on the main-sequence. Methods: We construct a rotational evolution model that is driven by observed rotational distributions of young stellar clusters. Fitting the free parameters in our model allows us to predict how wind mass loss rate depends on stellar mass, radius, and rotation. We couple the results to the wind model developed in Paper I of this series to predict how wind properties evolve on the main-sequence. Results: We estimate that wind mass loss rate scales with stellar parameters as $\dot{M}_\star \propto R_\star^2 \Omega_\star^{1.33} M_\star^{-3.36}$. We estimate that at young ages, the solar wind likely had a mass loss rate that is an order of magnitude higher than that of the current solar wind. This leads to the wind having a higher density at younger ages; however, the magnitude of this change depends strongly on how we scale wind temperature. Due to the spread in rotation rates, young stars show a large range of wind properties at a given age. This spread in wind properties disappears as the stars age. Conclusions: There is a large uncertainty in our knowledge of the evolution of stellar winds on the main-sequence, due both to our lack of knowledge of stellar winds and the large spread in rotation rates at young ages. Given the sensitivity of planetary atmospheres to stellar wind and radiation conditions, these uncertainties can be significant for our understanding of the evolution of planetary environments.Comment: 26 pages, 14 figures, 2 tables, to be published in A&

Stellar Winds on the Main-Sequence I: Wind Model

Aims: We develop a method for estimating the properties of stellar winds for low-mass main-sequence stars between masses of 0.4 and 1.1 solar masses at a range of distances from the star. Methods: We use 1D thermal pressure driven hydrodynamic wind models run using the Versatile Advection Code. Using in situ measurements of the solar wind, we produce models for the slow and fast components of the solar wind. We consider two radically different methods for scaling the base temperature of the wind to other stars: in Model A, we assume that wind temperatures are fundamentally linked to coronal temperatures, and in Model B, we assume that the sound speed at the base of the wind is a fixed fraction of the escape velocity. In Paper II of this series, we use observationally constrained rotational evolution models to derive wind mass loss rates. Results: Our model for the solar wind provides an excellent description of the real solar wind far from the solar surface, but is unrealistic within the solar corona. We run a grid of 1200 wind models to derive relations for the wind properties as a function of stellar mass, radius, and wind temperature. Using these results, we explore how wind properties depend on stellar mass and rotation. Conclusions: Based on our two assumptions about the scaling of the wind temperature, we argue that there is still significant uncertainty in how these properties should be determined. Resolution of this uncertainty will probably require both the application of solar wind physics to other stars and detailed observational constraints on the properties of stellar winds. In the final section of this paper, we give step by step instructions for how to apply our results to calculate the stellar wind conditions far from the stellar surface.Comment: 24 pages, 13 figures, 2 tables, Accepted for publication in A&

Stellar activity and planetary atmosphere evolution in tight binary star systems

Context. In tight binary star systems, tidal interactions can significantly influence the rotational and orbital evolution of both stars, and therefore their activity evolution. This can have strong effects on the atmospheric evolution of planets that are orbiting the two stars. Aims. In this paper, we aim to study the evolution of stellar rotation and of X-ray and ultraviolet (XUV) radiation in tight binary systems consisting of two solar mass stars and use our results to study planetary atmosphere evolution in the habitable zones of these systems. Methods. We have applied a rotation model developed for single stars to binary systems, taking into account the effects of tidal interactions on the rotational and orbital evolution of both stars. We used empirical rotation-activity relations to predict XUV evolution tracks for the stars, which we used to model hydrodynamic escape of hydrogen dominated atmospheres. Results. When significant, tidal interactions increase the total amount of XUV energy emitted, and in the most extreme cases by up to factor of $\sim$50. We find that in the systems that we study, habitable zone planets with masses of 1~M$_\oplus$ can lose huge hydrogen atmospheres due to the extended high levels of XUV emission, and the time that is needed to lose these atmospheres depends on the binary orbital separation.For some orbital separations, and when the stars are born as rapid rotators, it is also possible for tidal interactions to protect atmospheres from erosion by quickly spinning down the stars. For very small orbital separations, the loss of orbital angular momentum by stellar winds causes the two stars to merge. We suggest that the merging of the two stars could cause previously frozen planets to become habitable due to the habitable zone boundaries moving outwards.Comment: Accepted for publication by A&

The origin of organic emission in NGC 2071

Context: The physical origin behind organic emission in embedded low-mass star formation has been fiercely debated in the last two decades. A multitude of scenarios have been proposed, from a hot corino to PDRs on cavity walls to shock excitation. Aims: The aim of this paper is to determine the location and the corresponding physical conditions of the gas responsible for organics emission lines. The outflows around the small protocluster NGC 2071 are an ideal testbed to differentiate between various scenarios. Methods: Using Herschel-HIFI and the SMA, observations of CH3OH, H2CO and CH3CN emission lines over a wide range of excitation energies were obtained. Comparisons to a grid of radiative transfer models provide constraints on the physical conditions. Comparison to H2O line shape is able to trace gas-phase synthesis versus a sputtered origin. Results: Emission of organics originates in three spots: the continuum sources IRS 1 ('B') and IRS 3 ('A') as well as a outflow position ('F'). Densities are above 10$^7$ cm$^{-3}$ and temperatures between 100 to 200 K. CH3OH emission observed with HIFI originates in all three regions and cannot be associated with a single region. Very little organic emission originates outside of these regions. Conclusions: Although the three regions are small (<1,500 AU), gas-phase organics likely originate from sputtering of ices due to outflow activity. The derived high densities (>10$^7$ cm$^{-3}$) are likely a requirement for organic molecules to survive from being destroyed by shock products. The lack of spatially extended emission confirms that organic molecules cannot (re)form through gas-phase synthesis, as opposed to H2O, which shows strong line wing emission. The lack of CH3CN emission at 'F' is evidence for a different history of ice processing due to the absence of a protostar at that location and recent ice mantle evaporation.Comment: 10 Pages, 8 figures, Accepted for Astronomy and Astrophysic